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Abstract 

Recent theoretical work by Hall & Seddougui (1989) has shown that strongly non- 
linear, high wavenumber Gortler vortices developing within a boundary layer flow are 
susceptible to a secondary instability which takes the form of travelling waves confined 
to a thin region centred at the outer edge of the vortex. This work considered the 
case in which the secondary mode could be satisfactorily described by a linear stability 
theory and in the current paper our objective is to extend this investigation of Hall 
& Seddougui (1989) into the nonlinear regime. We find that at this stage not only 
does the secondary mode become nonlinear but it also interacts with itself so as to 
modify the governing equations for the primary Gortler vortex. In this case then, the 
vortex and the travelling wave drive each other and, indeed, the whole flow structure is 
described by an infinite set of coupled, nonlinear differential equations. We undertake 
a Stuart-Watson type of weakly nonlinear analysis of these equations and conclude, in 
particular, that on this basis there exist stable flow configurations in which the trav- 
elling mode is of finite amplitude. Implications of our findings for practical situations 
are discussed and it is shown that the theoretical conclusions drawn here are in good 
qualitative agreement with available experimental observations. 


x This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS 1-1 8605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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§1. Introduction. 


Our concern is with an asymptotic description of the three-dimensional break- 
down of steady, span wise periodic, large Gortler vortices via nonlinear travelling waves 
and is a natural extension of the work of Hall &; Seddougui (1989) who investigated 
the current problem in the context of breakdown through very small waves. 

Experiments by Dippes & Gortler (1972) and Aihara & Koyama (1981) suggest 
that the breakdown of vortices leads to a flow with wavy boundaries and which is 
periodic in time. Later, Kohama (1987) studied vortex instabilities in boundary layers 
over a laminar flow' control wing, and observed that a secondary instability of Gortler 
vortices appeared which was confined to a thin zone at the top of the extent of vor- 
tex activity. Additionally, he noted that the instability propagated with a speed that 
approached the free stream speed as the disturbance travelled downstream, and the 
onset of the unsteady flow perturbations was eventually followed by transition to tur- 
bulence. Peerliossaini & Wesfreid (1988) studied the fluid flow in a boundary layer 
in a concave section of a curved channel. Like Kohama, they found that when the 
secondary breakdown of Gortler vortices appeared it did so at the top of the vortices. 
Further dow'nstream they reported another instability, but which was confined to a 
zone close to the wall. 

Much of the early theoretical work concerned with Gortler vortices concentrated 
on the linear stability of external flows over concave walls. Notable contributions were 
made by, among others, Gortler (1940), Smith (1955), Hammerlin (1956) and Floryan 
& Saric (1979). However, Hall (1982a, b, 1983) demonstrated that much of this early 
work was fundamentally flawed for all the previous analyses had invoked the parallel 
flow approximation (which essentially assumes that the basic flow in which the vortices 
lie is independent of the streamwise co-ordinate and so neglects the effect of bound- 
ary layer growth). This approximation enables the linear instability equations to be 
expressed as ordinary differential equations but Hall illustrated that this assumption 
is unjustifiable except in the limit of small vortex wavelength, and indeed this is the 
explanation for the considerable inconsistencies in the results of the previous studies. 
Additionally, in the case of small vortex wavelength the Gortler instability may be 
described by an asymptotic structure which takes account of boundary layer growth 
in a fully rational manner and thence the parallel flow approximation is rendered su- 
perfluous in the only situation in which it has any relevance whatsoever. For larger 
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wavelengths the problem is fully non-parallel and the linear stability equations, which 
now take the form of a set of partial differential equations, have to be solved numeri- 
cally, see Hall (1983). This paper showed two significant features of this non-parallel 
flow problem, namely that the ideas of a unique neutral stability curve and of unique 
growth rates at a specified downstream location are inapplicable to the Gortler prob- 
lem because the location where a vortex commences to grow is dependent upon the 
position and the size and shape of the imposed disturbance. 

Questions concerning the development of nonlinear non-parallel vortices within 
growing boundary layers were answered by Hall (1988). This numerical investigation 
showed that as the nonlinear disturbance evolved, the perturbation energy became 
concentrated in the fundamental and mean flow correction; a conclusion consistent 
with the weakly nonlinear theory of Hall (19826) valid for small wavelength vortices. It 
is well known that Gortler vortices set up in an experiment conserve their wavelength 
as they move downstream. Since the boundary layer itself thickens it follows that 
the local nondimensional vortex wavenumber becomes large as the vortex develops. 
Thus the small wavelength limit in the external Gortler problem is appropriate to 
the ultimate development of any fixed wavelength vortex, and hence sufficiently far 
downstream in many flows the asymptotic work of Hall (1982a, 6) becomes applicable. 

As with all weakly nonlinear investigations, the results of Hall (19826) are valid 
only within the neighbourhood of the point where the considered perturbation is neu- 
trally stable. For vortices of wavenumber e” 1 , 0 < e <C 1, their development down- 
stream of the point of neutral stability is dictated by the solution of a pair of coupled 
nonlinear partial differential equations which adopt a simple asymptotic structure at 
large values of X, where eX (X= 0(1)) denotes the distance of the vortex down- 
stream of the neutral point. Formally, for large X , Hall & Lakin (1988) showed that 
this asymptotic structure could be used to deduce the flow configuration for fully 
nonlinear Gortler vortices, at which point the mean flow correction generated by the 
presence of the vortices becomes as large as the basic (undisturbed) flow itself. These 
fully nonlinear vortices, as described by Hall Sz Lakin (hereafter referred to as HL) 
are of the type whose stability to travelling wave disturbances will be considered here. 
For this reason, it is convenient to now briefly describe the main results of HL. 

Consider a vortex of wavelength 2ttc developing in a boundary layer flow at a 
Gortler number 0(e - ^), and suppose that the flow is neutrally stable r« i the down- 
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stream position x n . For x > x n , HL demonstrated that the flow field assumes the 
structure illustrated in Figure 1, where the vortex activity is restricted to region I and 
decays exponentially within thin shear layers IIa,b. Suppose that these are located at 
y i and j/2 * Then exterior to the shear layers, in zones IIIa,b, the vortex is annihilated 
and the mean flow satisfies the usual boundary layer equations. In the central part of 
the flow, region I, the mean flow is dictated by a solvability condition on the funda- 
mental component of the vortex and adjusts itself so as to make the fundamental and 
all harmonics finite within I. The vortex is determined by the mean flow equations, 
and the positions of the shear layers t/i,?/2 change as they move downstream: in fact 
they are given by a free boundary value problem derived from the boundary layer 
equations. Recall that this small wavelength asymptotic approach is relevant to many 
flows of practical importance for the reasons given previously. 

We shall study secondary instabilities of the flow in the shear layers IIa,b. Per- 
turbations which are spanwise periodic travelling waves are imposed on the flow within 
these zones and their development monitored. These disturbances are 7 t/ 2 radians out 
of phase with the fundamental in the spanwise direction, so that any secondary insta- 
bilities which occur produce locally wavy vortex boundaries in IIa,b. At the outset 
it is not clear that the shear zones should be susceptible to the wavy modes but Hall 
& Seddougui (1989) (hereafter referred to as HS) presented an argument as to why 
this should be the case. Basically, they considered a model equation which governs 
the nonlinear growth of time-dependent Gortler vortices in curved channel flows. The 
instability of a steady, streamwise independent vortex to travelling vortex-like pertur- 
bations was investigated and it was demonstrated that if the travelling perturbations 
were in phase with the fundamental component of the steady vortex then the flow is 
always stable whereas if the disturbance were 7 t/ 2 out of phase then instability could 
result. Moreover, if the Gortler number was large, the disturbance was confined to 
the shear layers and decayed away from the centres of these layers. This analogy with 
the curved channel flow problem led HS to analyse disturbances in the boundary layer 
problem which were concentrated in the shear layers Ila, b. As the disturbances de- 
cayed exponentially away from the centre of these shear layers, HS deduced that unless 
exponentially small terms were to be matched the two shear layers could be treated 
independently. 

In HS the wavy perturbations were assumed to have an infinitesimally small am- 
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plitude. It was found that in this case these wavy disturbances are fixed by a pair 
of coupled, second order, linear ordinary differential equations which were solved nu- 
merically. Hall & Scddougui found only one neutrally stable travelling wave solution 
of these equations and this information was used to explain some of the experimental 
observations of wavy vortices; notably those of Kohama (1987) and Peerhossaini & 
Wesfreid (1988). However, later studies by Bassom &: Seddougui (1989) modified the 
conclusions of HS by demonstrating that infinitely many neutrally stable modes are 
possible and properties of these modes in a high-wavenumber limit were described 
asymptotically. 

Here we aim to increase the amplitude of the wavy disturbances to extend the 
ideas of HS to the nonlinear regime. Importantly, it is found that at the point at 
which the wavy mode becomes nonlinear, the steady vortex flow as determined by 
HL is affected by self-interactions of the -wavy mode. Additionally, at this stage 
all the harmonics of the wavy mode become sufficiently large so as to influence the 
fundamental at leading order. Consequently, we found that this nonlinear problem is 
governed by a fully coupled, infinite set of differential equations which necessarily have 
to be solved numerically. Rather than tackle this complete system, we concentrate on 
a weakly nonlinear limit which enables further analytical progress to be made. The 
results of this work are considered in the light of the available experimental evidence 
and conjectures made concerning the results which could be anticipated from a strongly 
nonlinear calculation. 

The procedure adopted in the remainder of the paper is as follows. In section 2 
we formulate the problem and develop the defining equations for the nonlinear wavy 
modes in secti on 3. We analyse the small amplitude limit of these equations in section 
4 and in section 5 we discuss our results, consider the full numerical solution of the 
problem and draw some brief conclusions. 


§2. Formulation of the problem. 

The flow considered here is that as described in HL. We study the flow of an 
incompressible, viscous fluid of density p and kinematic viscosity v over a wall of 
varying curvature - ~x (A”/ L). Here R is a typical radius of curvature, X denotes 
distance along the wall and L is a typical lengtlxscale along the wall. Much of the 
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following detail concerning the basic flow structure in the absence of the imposed 
travelling disturbances follows very closely that described in HL and HS and so we 
keep our description of this aspect as brief as possible consistent with the aim of 
retaining clarity. 

The Reynolds number Re for the flow is defined by 


Re = 


U 0 L 

> 

V 


( 2 . 1 ) 


where Uq is a typical fluid velocity. The curvature parameter 6 is defined by 


<5 = 


L 

R' 


( 2 . 2 ) 


and our concern is with the case in which Re — > oo whilst the Gortler number G, 
given by 

G = 2 Re* 6, (2.3) 


maintains a constant 0(1) value. We denote time by T and co-ordinates along the wall, 
normal to the wall and in the spanwise direction by X, Y and Z respectively. If the 
corresponding velocity vector is (17, V , W) then we choose nondimensional co-ordinates 
(x,y,z) and velocities (u y v,w) given by 

(x,y,z) = L~ 1 (x,Re’Y,Re*z'j , (u,v,w) = Uq 1 (u, Re±V, Re^w) . (2.4a, 6) 


In the work considered here we restrict ourselves to examining flows for which 
u — ► 1 as y — > oo, although there is no formal difficulty in extending our treatment 
to account for more general flows. We write the pressure 


P = 



(2.4c) 


and substitute the expressions (2.4) into the governing continuity and Navier-Stokes 
equations. Defining the nondimensional time variable by t, where t = UqL~^T, we 
find that when terms of order Re~^ are neglected in the governing equations, these 


equations assume the forms 


du dv dw 

dx ^ dy dz 


du du du du d 2 u d 2 u 

dt + U dx + V dy + W dz dy 2 dz 2 ’ 


(2.5a) 

(2.56) 
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(2.5c) 


dv dv dv dv 1 o dp d^v d^v 

<9x dy dz 2 c?y dy 2 <9^ 

c?iu dw dw dw dp O^w d^w 

~m +u a7 + v 'd^ + w ^ = ~Yz + W + 'a^' 


(2.5 d) 


The motivation for the work contained in HL was to obtain a solution of (2.5) 
which satisfies 


u — v = w = 0 on y = 0, u ► 1 as y — > co. (2.6) 

and which is an asymptotic solution valid in the limit of small vortex wavelength. The 
Gortler number G expands according to 

G = Got ^ + G\e ^ + , , . , (2-7) 

where e~ l (0 < € <C 1) is the wavenumber of the vortices. In the core part of the 
boundary layer (region 7), the leading order 2 -independent part of the downstream 
velocity component u{x,y) satisfies 

G ° XU dy = lj (2.8) 

and this flow is driven by the finite-sized steady vortices. These vortices produce a 
mean flow correction term which is as large as the basic flow which would exist in the 
absence of the vortices and so the total mean flow field is altered at leading order by 
the presence of the vortices. The calculation for the vortices within the core region 
(see Figure 1) reveals that the vortices are trapped between yi and y2 (the values of 
which are found in the course of the computation) and decay to zero within the thin 
shear layers 77a, 5. 

The work given by Davey et. al. (1968) showed that Taylor vortex flow is unstable 
to perturbations differing in phase from the fundamental component of the steady 
vortex flow by 90°. In HS, the authors sought a secondary instability that produced 
locally wavy vortex boundaries in the shear layers 77a, b; akin to the wavy surfaces 
predicted by Davey et . al. in the case of the secondary instability of their Taylor vortex 
flow. In keeping with this aim, HS examined solutions for out of phase perturbations 
proportional to 

( i f x i£lt\ 

E = exp - ] K ( x)dx J , (2.9a) 
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where the wavenumber K was expanded as 

K = K 0 + e>K i + ..., (2.96) 

and fl is the constant frequency. The wavy vortex mode was assumed to be a short 
wavelength, high frequency mode which was necessary in order to ensure that this 
secondary instability was to remain trapped within Ila or 176, an assumption which 
is consistent with the available experimental observations. The scales in (2.9) are 
chosen, from Hall (19826), so that ~ ^ and f;- + it£- = 0 in the shear layers. The 
latter scaling ensures that the waves travel with the speed of the fluid in the shear 
layer, as observed in the work of Kohama (1987). 

In HS the secondary instabilities were supposed to be of infinitesimal amplitude 
so that the steady vortex flow was unaffected by the presence of these wavy modes. 
Our present aim is to allow this travelling vortex to be of sufficiently large amplitude 
so that self-interactions of this mode do modify the governing equations for the steady 
vortex. We then have the case of a fully coupled problem in which the primary vortex 
and the wavy mode must be determined simultaneously; in contrast to the calculation 
presented in HS in which the wavy mode was found only after the steady vortex had 
been fully determined. In the following section we develop the asymptotic expansions 
for both the vortex and the nonlinear wavy mode under the assumption that the latter 
is confined to within the shear layer Ila, an analysis which may be trivially modified 
to account for modes within the other shear layer if so desired. 


§3. The asymptotic structure of the wavy modes. 


In HL it was demonstrated that the layers IIa,b are of thicknesses 0(e$) and so 
in Ila we define the co-ordinate 

(y - V2 ) 


£ = 


3 _ 5 

es 


(3.1) 


where y 2 (z) is the location of the zone Ila. Hence, in this layer we replace d/dx and 
d/dy by 

d y 2 d Id 

t — and — — - 

dx di ef di 

respectively. The vortex activity is wholly within the core region I between Ila and 
176, and the algebraically decaying vortices in I are reduced to zero exponentially 
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within the shear layers. In the analysis described below we find that the leading order 
mean flow terms across the whole of the boundary layer are unchanged by the presence 
of the nonlinear wavy mode terms. Additionally, the steady vortex terms within the 
core region I are unaffected by the wavy mode, but this latter disturbance does alter 
the vortex quantities within the shear layers. 

We consider the amplitude of the wavy mode to be such that the equations given 
in HS to determine the steady vortex within the shear layers are modified due to 
interactions of the wavy disturbances. It is found that when the wavy mode is as large 
as the steady vortex (= 0(e * ), by HL) this change in the character of the governing 
equations occurs and we are led to consider total velocity and pressure fields within 
the shear layer Ila which take the forms 


u =uq + € * Si + eU 2 + f 3 83 + . . . + €* [cu 0 '0 + e 3 ?CU 0cl + ...} 
+ e 


E (CU0cn + SuO.n) E (C«lcn + Sll\, n ) E" + . 

n— 1 n— 1 


. + c.c. 


+ 


u 0 nE n 4- . . . + c.c. 

n —1 


+ * • • j 


(3.2a) 


v =ao + + • • • + c 3 [cVo c o + e t CVq c \ + . . .] 

+ 


E (C^Ocn + Svq ,n)E n + C 3 X] (Cv l en + Sv\ sn ) E n + 

n =1 n — 1 


+ C.C. 


+ 


V On E 71 + . . . + c.c. 

n — 1 


+ « - - j 


(3.26) 


w 


=e 3 5TFo a o + € 3 £Wq»i + . . . 


+ e‘ 


+ e~ 


E (C^Ocn + Swo sn )E n + e 3 E (C w lcn + Swi, n )E n 4- . 

n =1 n —1 


+ C.C. 


E WOnE” + ... + C.C. 

n —1 


+ 


(3.2c) 
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P =e' 


' 4 Po + 


J Pi + 


+ e 3 CPq c q + e 3 CPod + 


+ e' 


YL ( CpOcn + SpQsn) E n + (Cpi en + Spi, n ) E* + . . . + C.C. 


1=1 


i=l 


+ 6 » 


y! P On E n + • • - + C.C. 
1 


+ .... 


(3.2 d) 

In (3.2) C and S denote cos(z/e) and sin(z/e) respectively, all unknown coeffi- 
cients are functions of x and f and c.c. denotes complex conjugate. The expansions 
are more complicated versions of those presented in HS, and we recall that the vari- 
ation E defined by (2.9) dictates the streamwise development of those parts of the 
total velocity and pressure fields corresponding to the wavy mode. As is expected, 
when the amplitude of the wavy mode is decreased from that value implied in (3.2), 
we retrieve the expansions of HS. We remark that higher harmonics of the spanwise 
variation lead to terms of smaller sizes than those given in (3.2), see HL. These terms 
do not enter the determination of the primary vortex or the wavy mode and so have 
been omitted in (3.2). However, we do observe that all harmonics of the spatial and 
temporal variation E must be accounted for since they all play nontrivial roles in the 
subsequent analysis. 

We find the unknowns in (3.2) by substituting these forms into (2.5) and equating 
coefficients of specific terms. The leading order mean flow quantities So and Si are 
found to be precisely as in HL so that 


a + 2y 2 
Gox ’ 


V Gox(a + 2 y 2 ) 


(3.3) 


Here a(x) is an arbitrary function which may be determined by considering the so- 
lutions in region J, see HL, In HL, the fundamental vortex term Vo c o was found by 
deriving a solvability criterion on higher order vortex terms within the shear layer 
Ila , With the large amplitude wavy disturbance imposed, this solvability condition 
assumes the form 


o d 2 V 0c 0 £Vo c o /— — 7 — — — _ 

3 - 2 1 T" n *■ V ^0X( a + 2j/2)-^7-^0e0 - 

0^ a + 2 y 2 oti 

oo oo 

E ('i'on Oo Jn - Goxwowo.n)* + c.c.) - GoXMoXl 


.=1 


»=1 


V0c n -^r + c.c.), 


(3.4a) 
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where 


du 3 t , \ 

= fix) 

di K 


X(a + 2 y 2 ) 


3GoX‘ 


+ 


fa + 2?/o) ' 


t/2 

^OcO 


\/ G , ox(« + 2^2)^ 2\/ Gox(« + 22/2) 


+ 


s/(£ 

n =1 


Ocn , f|A 

a“ °' ,_ ar 


+ WOcnW5, n - W0,n«0 C n + C C • J I ^ 


(3.46) 

These equations are derived following procedures identical to those described in 
HS. Here b(x ) is a function arising from the solution in the core, f(x) is a function 
which can only be found by considering higher order terms (which proves unnecessary 
here) and an asterisk on a quantity denotes the complex conjugate of that quantity. 
The equations equivalent to (3.4) for the case of infinitesimally sized wavy vortices, 
as in HS, may be retrieved by removing the terms comprising products of wavy mode 
perturbation quantities in (3.4a, 6). 

Comparing coefficients of sin (z/e) in (2.5a — c) and of cos (z/e) in (2.5d) yields 

the following relations between the leading order wavy vortex terms. We obtain 

dvo, n 


dZ 


- men = 0, 


in (Kquq — £1) MOin + V0,n 


du 1 


= -«o. 


in (Kquq - £1) l>o, n + GoXW 0 «O.n = — v 0.n, 
in (K 0 Q 0 - fi) tnocn = ~P0.n ~ WQ cn . 
Now (3.56, c) are only consistent if Kquq = O and then we have 


WOcn = 


0VQ sr 


du\ 

~aT' uo,n = 

By an entirely similar procedure we can obtain 

dvo cn du\ 


VQ.n, 


PO.n = 


m, r 


UQcn 


-VQcn , 


POcn 



(3.5a) 


(3.56) 


(3.5c) 


(3.5d) 

9v 0 , n 

’ 

(3.5e) 

dVQen 

dt ’ 

(3.5/) 


Substituting into (3.4) and simplifying reveals that the steady fundamental vortex 
term T^ c q satisfies “ 


0 2 y OeO 

di 2 


<Vq c q 


x (a 4- 2y2) ; 
9X\/^0X 


6 _ 1 
3 3(a + 21 / 2 ) J 


t/3 

v 0c0 

6 


1 1 (k.i 2 + K.P) ) - f K^KL±M Vo .o 


n =1 


(3.6) 




dU r 


2 W Q n UQa n — V 7 Gox ( a + 2) q, ^Oen + c - c 
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Examining terms proportional to E n in (2.5 d) yields the governing equation for 
the mean flow correction term W o n '> 


PWOn 

dZ 2 

I ( v '°-° 


— in 


T/ - = . ^^1 

K\uq H — =■ 


5 2 uq,, 


3* 


2 ^0, n ’ 


u 0 

S^VocO 

a^ 2 


VEOn = 


1 v d f ^^O.(n-k) 0c (n - fc ) A 

J+5 E aj ”»•*— — )■ 

Jfc =— CO 


(3.7) 

To obtain the equation for the wavy vortex term uo,„ we need to consider the 
governing equations for the higher order terms in the scheme outlined in obtaining 
(3.5). Following this procedure we find from (2.56, c) that ui, n and i>i, n satisfy 

d 2 uo, n 


in (A'ofti + Kluo)uo, n + UQjn "^7 ^On^OcO 


di 


^ ( dU 0{n _ k) jxr 

+ Lj l v 0.k — WokUQ, 


jfc =- CO 


di 


c(n-fc)^ — u l»n U l.n 


du 1 

aT 


(3.8a) 


and 


d 2 vo, n , dpo„ 


in ( Kqui + ATiuo) vo. n 7^2 


~ VbcO^On + GoX^lWO.n 


VQck 1^0(n- k) — — GoX U O u l«n W l*n ■ 


it =— OO 


(3.86) 


Now using the results (3.5 e,f) we find that (3.8a, 6) are consistent for the unknowns 

(«l*n,n.n) ° nl y if 


a 2 up <n 

dZ 2 


+ £ 


1 x' (q + 2y2) r _ 6 

3(a + 2 t/ 2) 9x\/£oX 3 J 


VQj 


2m 

~3~ 


Kly/a + 2 1/2 „ 1 

+ . ^ j U 0.n = 


^Ojn^OcO 


>/Gox a + ^2/2^ 


+ (ko cn | 2 + ko, n | 2 ) - 7 V G 0 x{a + 2y 2 )v 0 . n ~ |lU 0 n U 0c o 

3 n=l 3 

— — ^ \/G'0X( C1 + 2y2) u 0.k H 2t>0ck ^0(n- k)^ • 

jfc — - 00 


2.:. 


(3.9a) 


Using an identical method to determine a solvability condition for the quantities 
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u\ cn and v\ cn , we also find that 


d 2 V0cn 

oe 

2m 

T” 




1 X_ (a + 2my _ b_ 

3 (a + 2y 2 ) 9Xn/^0X 3 


VQc 


Kl\/a + 2 j/2 

+ 


y/GoX 


a + 2j/2 J 


^Ocn 


^OenVr 2 


OcO 


+ — t— X] (K’Ocnl 2 + N,„| 2 ) - {rV Gox{a + 2 y 2 )v 0cn - ~V Gqx(« + 2 y 2 )‘^r 


X=1 


dZ 


- - S yV Gox(a + 2y 2 )vo 

k — — oo 


^^O(n-k) T J 7 | 

ct ^ 2U0 f fcW0(n-t) J • 


(3.96) 

Filially, to close the system, we need to find the equations which determine the 
harmonic terms ?7o n . Considering the coefficients of E n in (2.5a, b), we find that 


and 


a 2 c/o„ 

di 2 


— in ( 


■ v h , dV ^ n 

IuRqU On H T— = 0, 

di 


ftui \ * dui * Id 

K iSo + -r 1 ) U 0n ~ ~^V 0n = - = — 

^0 7 v Gqx (° + 2y 2 ) 


(3.10a) 


(^OerVOcO) 


d 


(vo ek v 0c ( n _ k ) + uo.it no. (n-*)) 


2v / G 0 x(a + 2y 2 )^ L* = _ ( 

(3.106) 

We now have the complete specification of our problem to determine the nonlinear 
vortices. The steady vortex term Vo c o is given by (3.6), the mean flow Won by (3.7), the 
wavy mode quantities v$ en , i>o, n by (3.9) and the harmonic terms U o n by (3.10). This 
system comprises an eigenproblem which may be simplified by making the following 
transformations. We define the functions 


Ql(x) - 


1 


6 X (a + 2y 2 ) ] 


3 (a + 2y 2 ) 3 OxV^oX 

transform the independent co-ordinate £ to £1 by defining 

92 f 


\/Gox(a + 2y 2 ) 

g 2 (x) = , (3.11a, 6) 


ft = £ + 


write 


A'J = K] 


91 

fg 2 9l\/Gox 

I 

gi (a + 2 \y 2 ) 2 


(3.11c) 


(3. lid) 
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and finally define nondimensional quantities according to 


=(-<7l)*fl> Vb c o = \/6(-0 i)»A, vo cn - (~gi)‘Bo cn , vo. n = (-gi) 3 Bo. n , 
t"t7 / 9 ^0n (-9lY T „ o n 

'V On = {-gi ) 3 Con, — = D On, “ = 


dg 


92 


(“5l)(« + 2y2)’ 


and 


T y \/a + 2j/2 , 2. r .f 

Al = ^ r ( " 9l) Kl 


(3.12) 


Inspection of the resulting equations reveals that necessarily we must have 


#0c(2;-l) = Bq 8 (2j ) = Co(2;) = -Do(2j-1) = 0, 


for all integers and then if we define 


Bo c ( 2 j) = B 2 j, Bq, ( 2 j _ l) = Q)(2;- 1) = Cj , D 0 ( 2 j) = Dj, (3.13a) 

note that by the conjugacy relationships implied by (3.2) we have 


B-j^B', D_ j= D;, C } = Ci_j, 


(3.13b) 


and define the quantities B 0 = Dq = 0 for completeness, we obtain the system 

d 2 A 


dg* 


— g A =A J + 




(3.14a) 


+ — 7=5^ \ 2 B 2 i-xCx-i - 3B2jD_j + c.c.] , 

3v/e , =1 


- D %~ 1 ~ gB 2 j-\ - 2 Yr(2j - + Ki)B 2j -i = B 2j -iA 2 


dg* 


+ 


B2j - 1 


E | b,i 


2\/6 


- V B 21 - £(2i)(ftt| + = B 2j A 2 

dr ) z 3 


3 CjA ~\ ^ (3B2r-\Dj- r + 2B2 T Cj . T ) , 

r 00 

(3.14b) 


+ 


^-(e |B.| 2 ) -AD^-l E (3B 2 r fl j . r -2B 2 r _ 1 C,_, +1 ), 

r— 1 r — — 00 


(3.14c) 
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(3.14c?) 


and 


d l Cj 


V - i(2j - !)(&, + Kl)C, = f - ^- 1 0) 


d 2 Py 


dry 


j- — 2ij(Clr) + I<\)Dj = 


1 d 4 


3 drp L 2 _ 


1 °° 

“ X] {B 2r B2j-2r + -B2r-l-02j_2r4-l) + y/sBzjA 


(3.14e) 


As yet we have not addressed the problem of specifying the boundary conditions 
pertinent to the system (3.14). Firstly, we anticipate that the wavy vortex be con- 
strained to lie completely within the shear layer I la and this requires that Bj , Cj and 
Dj decay exponentially as |?7 1 — ► oo. Secondly, the steady vortex term A must match 
with the core region solution for the strongly nonlinear vortex within I. This, as in 
HL, implies that A ~ (— rj)? as rj — ► — co. Finally, all the vortex activity is confined 
within the core I and the shear layer Ila so that in order that in Ilia the mean flow 
is governed by the usual boundary layer equations we impose the requirement A — > 0 
as rj — > oo; indeed by (3.14a) we have A oc Ai(rj) in this limit where Ai is the Airy 
function, see Abramowitz &: Stegun (1964). 

To conclude this section, we note that the nonlinear wavy modes within the shear 
layer Ila are determined by the functions Bj , Cj , Dj which are defined in (3,12) and 

(3.13) . These terms, together with the steady vortex quantity A, satisfy equations 

(3.14) with associated boundary conditions 


A oc Ai{rj) as rj — > oo, (3.15a) 

A ~ (— rj)^ as rj — ► — oo, (3.155) 

and 

B j ,Cj,D j — > 0 as |t ?| — > oo. (3.15c) 

In general, the solution of this system is highly complicated numerical task and 
in the following section we address this problem in a weakly nonlinear limit which 
permits additional analytical progress. 
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§4. Weakly nonlinear theory. 


We recall that in HS the infinitesimally small wavy mode problem was governed 
by the system of equations (3.14) for the functions A, B\ and C\ when all the nonlin- 
ear terms (except those involving the function A) are omitted. HS used a numerical 
scheme to search for eigenvalues of this linearised problem but located only one such 
pair. However, HS suggested that further eigenvalue pairs with K\ and Cl real (and so 
corresponding to neutrally stable wavy modes) may be possible and Bassom & Sed- 
dougui (1989) found eight such pairs in all. Further, these latter authors also demon- 
strated that, plausibly, infinitely many pairs exist and they derived the asymptotic 
structure of the neutrally stable modes in a high-wavenumber, high-frequency limit. 
The eight eigenpairs identified by Bassom Seddougui (hereafter BS) are summarised 
in Table 1. 

Our analysis of the full equations (3.14) is concerned with the problem in which 
the wavy mode is of amplitude 0(h), ft < 1, relative to the scalings implied in (3.2). 
In this case, in keeping with the usual weakly nonlinear approach, we anticipate from 
(3.14) that the perturbation quantities are changed by 0(h 2 ) from their linear values, 
and we write 

A =Aq + h 2 A\ + . . . , B\ = HB\q + h 2 B\\ + . . . , Z ?2 = h 2 B20 + . . . , 

C\ =hC\Q + h 3 Cn 4- . . . , C 2 = h 2 C20 + . . . , D\ = h 2 D iq + . . . , (4.1) 

Ci = (h) L + h 2 n p + ..., K\ — (ki) L + h 2 K p . 

The remaining terms in (3.14) are sufficiently small so as to become negligible in the 
ensuing analysis. 

Substituting (4.1) in (3.14) and (3.15) gives, at leading orders in h, the linear 
eigenvalue problem described by HS. In this paper, together with BS, figures of the 
linear eigenfunctions B 10 and C 10 appropriate to some of the eigenpairs in Table 1 are 
shown. We notice from (3.14a), (3.15 a,b) that the steady vortex profile Ao satisfies 

a ,3 

= -v 

Ao — > \Ai(r}) as 77 — > co, Ao ~ (— 77)^ as 77 ► —00. 

Further, the results of Hastings & McLeod (1980) demonstrate that A = \/2 and then 
the function Ao may be evaluated by using the asymptotic condition as 77 — > 00 and 
integrating in the direction of decreasing 77 using a Runge-Kutta method, see HS. 
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At next order in h when (4.1) is inserted in (3.14), we obtain 

( J - 2 - T l) M - SAqAi = -A 0 |i?lo| 2 + —7= (-BioCio + C.C.) , 


where 


drj A 

fBn 

dr) 2 

d 2 C n 


3\/6 

+ ai(v)Bn + a 2 (v)Cn = Sl p £i(r}) + K p £ 2 (rj) + £3(^7), 

\ZqAqC\{t}) 


^2 — ^ G 3(^)-Bll + a 4 ( r ?)Cll = ftp 


ir)C 10 + 


+ K P 


iC W + 


V^AoC 2 (r}) 


+ ^£3(0) + CiM, 


(, , 2ih L ~] 2i(K 1 ) L a2 2Ve A 

ai = -7? I 1 + -y- J A 0 , a 2 = ~Y~ a 0 , 


(4.2a) 

(4.26) 

(4.2c) 


(4.3a, 6) 


a 3 


iy/c 


(fl L 71 + (Ki)l) Ao, 04 = -i (n L r) + {Ki)l) + 2Aq. (4.3 c,d) 


Here, we find that . . . , £4 are defined by 

2 i 2 i 

Cl = —77-010, 02 = —010, 


(4.4a, 6) 


£3 s 2 AoAxBio + glol ^ 1()|2 - - < ^-C 10 A 1 - B{ 0 D W - 2 -B 2Q C* w , (4.4c) 

£ 4 = - Vsa x aIb v 0 - 2A0A1C10 + l -~ (n L v + {ki)J AiBw - Z?10 ^ 1 °l 2 

- ^ (B 10 C{ 0 + 2B[ 0 C 10 ) - i (n LV + (Ki) ) B' w B 20 

v/6, „ , £>io^o-Oio\/6 

y-4o#2oCio H 2 ’ 

(4.4c?) 

where, from HS, we know that the functions 0 10 and Cio satisfy the homogeneous 
forms of equations (4.26, c). Additionally, from (3.14) we obtain 

— T)B 2 0 - %■ (&lV + (a"i) ) B 2 o - AqB 2 o = “^BioCio - VgDioAo, (4.5a) 
dr}* 3 l 3 


and 


d 2 D 


10 


1 d 2 


j-5- + (*l) )D W =- j S j)2 


j(Bio) 2 


+ V6i?20^o] ■ (4.; 
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To ensure the disturbance is confined within the shear layer we demand Ai, D n, 
C 11, 13 20 and D io » 0 as |r?| > oo. 

The homogeneous forms of (4.26, c) are solved in HS. Then, as is the normal 
procedure within a weakly nonlinear analysis, (4.26, c) only have a suitable solution 
if a certain compatibility condition is satisfied. This condition leads to the speci- 
fication of the correction terms K p and Q p within the wavenumber and frequency 
expansions detailed in (4.1). To derive our compatibility requirement we first need 
to consider the system adjoint to the homogeneous forms of (4.26, c). The adjoint 
functions {Gi(r/), G 2 (r])) satisfy the coupled equations 


d 2 G i 
di] 2 

d 2 G 2 

d 77 2 


+ ai(v)Gl + a3(r))G2 = 0 , 


+ (12{v)Gi + a^(r])G 2 = 0, 


(4.6a) 


(4.66) 


with associated boundary conditions of G i, G 2 — * 0 as |7/| — » oo. Multiplying (4.26) 
by Gi(r}), (4.2 c) by G 2 (ry), adding the resulting equations and integrating by parts, 
reveals that for a satisfactory solution of (4.26, c) to exist, we require 


J (GiNi + G 2 N 2 ) dr} = 0, (4.7) 


where N\ and N 2 are the right hand sides of (4.26, c) respectively. 

On substituting (4.4) into (4.7) we eventually obtain a relationship between the 
wavenumber correction K p and the frequency correction of the type 


1 + zifl p + z 2 K p = 0 , 


(4.8) 


where z 1, z 2 are complex numbers expressible as quotients of integrals. In particular, 
z\ = J1/J3 and z 2 — J2/J3 where 


h 


-r 


C\G\ + j G 2 


dr}, 


(4.9a) 


and 


h 


= /” [rsoi + (ic lc + g 2 


dr h 


h 


-r 


£ 3 Gl+ (v^ + A ) G2 


dr}. 


(4.96) 


(4.9c) 
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Before presenting results of the computations, we observe that we may easily 
generalise our analysis to allow the wavy perturbation to evolve on a slow lengthscale 
so that the growth or decay of a non-neutral disturbance may be monitored. Formally, 
if we consider the neighbourhood of a point x — x n (at which an infinitesimal wavy 
vortex is neutrally stable) and allow a constant frequency disturbance to develop on 
the lengthscale X , where X = e~ >h~ 2 {x - x n ), h < 1 , then the result of repeating 
the previous analysis for a disturbance of amplitude hA(X ) relative to the scalings in 
( 3 . 2 ) is the evolution equation 


dA 

lx 


Vl p A - — \A\ 2 A. 

Z2 Z2 


( 4 . 10 ) 


We may then easily obtain the equation for the development of the wavy disturbance 
amplitude 

A (lil 2 ) = 2 Re ( Al) a, |i| 2 + 2 Re ( = 1 ) |i| 4 . ( 4 . 11 ) 

dX ' ' ' Z2 ' 'Z2' 

To effect the necessary calculations we used the method described by HS to ob- 
tain the linear eigenfunctions D\q and C\q. To fix a normalisation for the problem, we 
elected to scale these eigenfunctions such that max{|Z?io| : |q| < 00} = 1. The tech- 
nique of HS is easily adapted to solve the adjoint system ( 4 . 6 ) and essentially consists 
of the following. Equations ( 4 . 6 ) were written as a set of four first order differential 
equations, which was solved by using a fourth order Runge-Kutta scheme. Since ( 4 . 6 ) 
is the system adjoint to the problem specified in HS, the two problems have the same 
eigenvalues and the integration procedures for ( 4 . 6 ) were started at 77 = —00 and at 
77 = 00 and both continued to 77 = 0. The detailed asymptotic conditions on G\ and 
C?2 for large [77 1 follow in a similar manner to those given in HS and so yield the ini- 
tial conditions for the integrations. The two integration procedures from ±00 were 
matched across 77 = 0, and hence the adjoint functions G\, G2 determined. 

In order to obtain the integrands in ( 4 . 9 ) we needed to find the functions A\, B20 
and Dio in addition to those already computed. The steady vortex quantity A\, which 

satisfies jyli| > 0 as |?7 1 > 00 was evaluated using a finite differencing procedure on 

(4.2aj. We observe from ( 4 . 5 ) that the remaining wavy vortex terms Dio an d B 20 are 
governed by coupled equations and these functions were found also by employing finite 
difference methods. Once these flow quantities were all determined, the necessary 
integrations to evaluate ( 4 . 9 ) were effected by use of trapezoidal integrations and 
the results so obtained were checked to ensure that they were independent of thv. 
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step lengths used within the Runge-Kutta and finite differencing procedures and also 
independent of the (large) values of |ry| at which the various asymptotic behaviours as 
Jt 7 | — ► oo were assumed to be valid. 

The results of these computations and their implications are considered in the 
following section. 


§5. Results and discussion. 


The numerical scheme outlined above was used to determine the coefficients in 
the amplitude equation ( 4 . 11 ). This equation demonstrates that the weak nonlinearity 
of the problem allows the existence of a threshold equilibrium amplitude A c given by 


|i. 


2 


Re (-izi z 2 | 0 

Re(-^ 1 ) ) ' 


(5.1) 


The relevant values of the quantities contained within (5.1) are summarised in Table 
2 , where we list those corresponding to the first five linearly neutrally stable modes 
given in Table 1. Eigenfunctions for the problem (4.2)-(4.5) in the cases tl L = 0.372 
and Ct L — 0.659 are shown in Figures 2 and 3 respectively. 

Inspection of the complex-valued eigenfunctions reveals that the imaginary parts 
are pictorially similar to the real parts and so we elect to concentrate upon the latter 
in presenting the eigenfunctions of the problem. Observation of the eigensolutions for 
the first five linearly neutrally stable modes shows that the eigenfunctions G j, C? 2 , ^ 20 > 
Dio have the same trends as the linear eigensolutions B\o, C 10 described by Bassom 
& Seddougui (1989). Thus, the eigenfunctions are oscillatory in nature and centred 
about a negative 77 position. Moreover, this position is close to r] = We 

also note that the number of oscillations of the complex eigenfunctions increases for 
each eigenmode with the eigenfunctions corresponding to the choice ((K i ) L , ftx.) = 
(0.690,0.372) having the least number of oscillations. 

Since the coefficient of the nonlinear term in (4.11) is negative for each of the 
cases considered in Table 2 , interest is concentrated on supercritical disturbances when 
ftp < 0 . Then, on the basis of linear theory, the wavy modes of Table 2 are unstable, 
but the nonlinearity in ( 4 . 11 ) ensures that the equilibrium amplitude A e is stable in 
the sense that for A < A t the disturbance will tend to grow whereas for A > A t it 
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will decay until A = A e . Of course, here we have only discussed the results for the 
five modes of lowest frequencies and from the work of Dassom &: Seddougui (1989) we 
know that it is probable that there are infinitely many more modes. However, it is 
the case that in the development of the flow over a curved wall the lower frequency 
vortices are likely to be the more important: an aspect which we consider in more 
detail later. 

The procedure described above sought a description of the evolution of a wavy 
mode in terms of the perturbation of the non-dimensional frequency away from the 
linear neutral values (fi) . We remark that in order to translate our findings for the 

L 

case of any particular curved wall of prescribed curvature x(:r), we may only determine 
the dimensional wavenumber K\ and frequency ft once the scaling quantities in (3.11) 
and (3.12) have been calculated. To do this requires an application of the computation 
described by HL for the particular x(z)- In addition to their account of this numerical 
calculation, HL investigated two asymptotic limits. Firstly, solutions were determined 
for their free-boundary problem for x close to the linear neutral position x = x. , and 
secondly, for x far downstream of this point. In the latter case it was assumed that 
the curvature x increases at least as quickly as for x 1. 

Rather than attempt to pursue the HL calculation for any specific boundary layer 
flows here, it was felt that it would be more instructive to draw general conclusions 
concerning nonlinear wavy vortices which would hold for large classes of flows. This is 
largely motivated by the fact that all the experimental papers relating to wavy vortices 
within boundary layer flows over curved walls contain insufficient detail to enable the 
dimensionless values of K\ and f l to be computed. Hence any results we may draw 
concerning specific flows could not be easily related to practical observations. 

We now turn to consider the behaviour of the frequency of the neutrally stable 
nonlinear modes in each of the asymptotic limits considered by HL. As x — > x , , 
where we recall the definition that x, is the location where the steady Gortler vortex 
is linearly neutrally stable, we know that the two shear layers merge. Then, following 
the notation of HS, if we denote by f l T and the (dimensional) frequencies of the 
neutral modes in the upper and lower shear layers respectively, we have from (3.12) 
that as x — > x, , — > Vl b . We now suppose that for large x, x ~ x M , M > j. This 

choice of curvature is again chosen to illustrate the type of behaviour we may observe 
in a practical situation rather than as representing any one particular boundary layer 
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flow. In this case, the behaviour of the functions a, yi, y2, b are as given by HL and 
by routine application of these forms, we may easily show, as in HS, that for x 1, 

n T ~ —mGqx 2 **- 1 , n B ~ (5.2) 

9 3 

Consequently, as x increases, the frequency of the neutral wavy mode confined to the 
upper shear laj^er increases whereas that of the mode in the lower layer tends to a 
constant value. Thus the modes can be categorised as being of high and low frequency 
respectively. 

We remark that the instability process detailed in section 3 for the upper layer 
is also perfectly possible within the lower layer //6, and if exponentially small terms 
may be neglected, the two shear layers are completely independent for the unstable 
travelling waves decay exponentially away from the centres of //a, 6. Theoretically at 
least, our analysis could be adapted to prove that an 0(1) disturbance in one shear 
layer provokes an exponentially small response in the other. 

In a practical situation, it is often found that a wavy mode of a predetermined 
frequency is imposed upon the flow and the development of the disturbance monitored 
as it moves downstream. Then we anticipate that the two shear layers will breakdown 
at different x -locations. As discussed by HS, since the downstream velocity component 
of the the basic state is larger within the upper layer //a, it is expected that this layer 
is the first to become unstable. Further, as the wavenumber K(x) relating to a fixed 
frequency perturbation will be different in the two shear layers, the modes in the layers 
move with different wavespeeds. Close to the linear neutral position, the shear layers 
meet and so if breakdown occurs close to this point then the structures within the two 
layers will be very similar. Far downstream, the layer Ila moves into the freestream 
and the downstream velocity component tends to the freestream speed. The lower layer 
sinks towards the bounding wall so the fluid velocity there approaches zero. Thence 
if the stationary vortices develop over a sufficiently lengthy interval before breakdown 
occurs, the mode in the upper layer has speed close to that of the freestream, whilst 
the one in 1/6 has a very small speed. 

We now consider the implications of our analysis for an experimental setting. 
We note from (5.2) that for the case of a mode of prescribed dimensional frequency ft 
imposed upon the flow upstream of the linear neutral point, as we move far downstream 
(for the case of x ~ ) the non-dimensional frequency ft for vortices confined to 
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the upper layer Ila tends towards zero. Our results obtained here demonstrate, as 
previously mentioned, the existence of a stable, threshold equilibrium amplitude, at 
least for the weakly nonlinear modes close (in parameter space) to one of the first 
five linearly neutral disturbances. As the perturbation propagates downstream, the 
non-dimensional frequency ft of modes trapped within the upper shear layer decreases 
and our analysis indicates that there are no unstable non-zero equilibrium amplitudes 
within this frequency regime. Hence, we would anticipate that the low-frequency (low 
f2) modes investigated in this paper are likely to be the most important in a practical 
experiment and that finite amplitude wavy vortices are likely to be observable in 
practice. 

In this paper we have chosen not to pursue an extensive solution of (3.14) away 
from the domain where the weakly nonlinear theory is applicable. In such a calcula- 
tion, all the unknowns in (3.14) are potentially important and would necessarily have 
to be accounted for in a fully rational way. Preliminary studies of the computation of 
the full system (3.14) and (3.15) suggest that this task is far from easy. It would be 
interesting to make a careful examination of these equations to determine the char- 
acteristics of strongly nonlinear wavy modes. However, here we have elected to deal 
only with the weakly nonlinear approach. The results of this analysis suggest that 
for numerical computations based on the full version of (3.14) it can be speculated 
that for each non-dimensional low frequency ft there is a multitude of possible vortex 
states, each stable and with a different wavelength. Equivalently, for a wavy vortex 
of prescribed streamwise wavenumber, there could conceivably be a family of possi- 
ble non-dimensional frequencies for this disturbance. A thorough description of the 
possible solutions of (3.14) is eagerly awaited. 

Papers by, among others, Bippes Sz Gortler (1972), Wortmann (1969) and Bippes 
(1978) describe some experimental investigations of the secondary instability of Gortler 
vortices. Recently, work by Kohama (1987) and Peerhossaini &: Wesfreid (1988) have 
given more details concerning this secondary breakdown than did the aforementioned 
authors. Kohama (1987) performed experiments on a NASA laminar flow wing. He 
considered the breakdown only in the upper parts of the vortices and gave measure- 
ments of the speeds of the wavy modes at various downstream locations. Peerhossaini 
&z Wesfreid (1988) examined a boundary layer over a concave section of a curved chan- 
nel. They reported that when the secondary instability first manifested itself, it did 
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so at the top of the steady Gortler vortices and that further downstream another sec- 
ondary mode was observed near the wall. Both instabilities had wavy boundaries and 
so, like HS, we conclude that these modes are formed as a result of the mechanism 
considered in this paper. This experimental work of Peerhossaini & Wesfreid (1988) 
lends credence to our conclusion that according to weakly nonlinear theory the finite 
amplitude wavy vortices are stable and hence expected to be obtainable in practice. 
Interestingly, Peerhossaini &: Wesfreid also commented that as the disturbance devel- 
oped downstream no coherent frequency of the mode was apparent. Possibly, in this 
case, the weakly nonlinear approach is no longer tenable and the full system (3.14) is 
appropriate; we then have the possibility of the existence of a whole family of modes 
with different frequencies as discussed above. This may explain the experimental ob- 
servation of Peerhossaini & Wesfreid (1988). 

To conclude, we observe that we have shown that according to our weakly non- 
linear theory, for non-dimensional frequencies close to the linear neutral values, finite 
amplitude, stable equilibrium wavy modes exist. Our results are consistent with the 
available experimental observations; however, the lack of detail given in the accounts of 
this practical work means that precise comparisons between those studies and our find- 
ings is rendered near impossible. Of interest would be the consideration of the strongly 
nonlinear system (3.14) although we ourselves are not planning to pursue this aspect 
further. However, we do feel that this investigation merits attention at some stage. 
We finally note that other instability processes could give rise to three-dimensionality 
and time-dependence in the latter stages of vortex development; these would most 
likely be Rayleigh instabilities associated with the spanwise locally inflexional velocity 
profiles, or Tollmien- Schlichting waves, both of which need to be thoroughly studied 
to enable a more complete picture of the transition process to be constructed. 
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mm 

0.690 

2.900 

4.156 

5.435 

7.53 



15.7 

■get 

0.372 

0.659 

0.742 

0.795 




1.60 


Table 1. Values of K\ and Cl for which infinitesimally small wavy 
vortex modes are neutrally stable, from Bassom & Seddougui (1989). 


m 


Kl 

— 

-2izi/z 2 

-2i/z 2 

EBBES1I 

0.6895 

0.3715 


(-0.187,0.345) 

3.89 

2.900 

0.6587 


(-1.019,-6.089 x 10- 2 ) 

2.033 

4.156 

0.7416 

(-5.64,6.70) 

(-1.044,-1.720) 

2.325 

5.435 


(-7.06,10.09) 

(-0.121,-3.997) 

7.652 

7.527 


(-4.57,11.86) 

(-1.335,-4.461) 

1.849 


Table 2. Quantities defined in (4.11) and (5.1) appropriate 
to the weakly nonlinear analysis of section 4. 


24 


ihhlillill UflllllU 































REFERENCES 


Abramowitz, M. & Stegun, I. A. 1964 A handbook of mathematical functions, 
Frankfurt; National Bureau of Standards. 

Aihara, Y. & Koyama, H. 1981 Secondary instability of Gortler vortices: formation 
of periodic three-dimensional coherent structures. Trans. Japan Soc. Aero. Space 
Sci., 24, 78-94. 


Bassom, A.P. k Seddougui, S.O. 1989 A note concerning the onset of three di- 
mensionality and time dependence in Gortler vortices. ICASE Rep. 89-77. Submitted 
to J. Fluid Mech. 

Bippes, H. 1978 Experimental study of laminar turbulent transition of a concave 
wall in parallel flow. NASA TM 7 5243. 

Bippes, H. k Gortler, H. 1972 Dreidimensionales storungen in der grenzscliicht 
an einer konkaven wand. Aera. Mech., 14, 251-267. 

Davey, A., DiPrima, R.C. k Stuart, J.T. 1968 On the instability of Taylor vortices. 
J. Fluid Mech., 31, 17-52. 

Floryan, J.M. k Saric, W.S. 1979 Stability of Gortler vortices in boundary layers. 
AIAA paper 79-1497. 

Gortler, H. 1940 Uber eine dreidimensionale instability laminare Grenzscliubten 
an Konkaven Wanden. NACA TM 1357. 

Hall, P. 1982a Taylor-Gortler vortices in fully developed or boundary layer flows: 
linear theory. J. Fluid Mech., 124, 475-494. 

Hall, P. 1982b On the nonlinear evolution of Gortler vortices in non-parallel 
boundary layers. J. Inst. Maths. Applies., 29, 173-196. 


25 


Hall, P. 1983 The linear development of Gortler vortices in growing boundary 
layers. J. Fluid Mech 130, 41-58. 

Hall, P. 1988 The nonlinear development of Gortler vortices in growing boundary 
layers. J. Fluid Mech., 193, 247-266. 

Hall, P. Lakin, W.D, 1988 The fully nonlinear development of Gortler vortices 
in growing boundary layers. Proc. R. Soc. Lond. A , 415, 421-444. 

Hall, P. Sz Seddougui, S. 1989 On the onset of three-dimensionality and time- 
dependence in Gortler vortices. J. Fluid Mech., 204, 405-420. 

Hammerlin, G. 1956 Zur Theorie der dreidimensionalen Instabilitat laminar Gren- 
schichten. Z. Angew. Math. Phys ., 1 , 156-167. 

Hastings, S.P. & McLeod, J.B. 1980 A boundary value problem associated with 
the Second Painleve Transcendent and the Korteweg-de Vries equation. Arch, ration 
Mech. Anal., 73, 31-51. 

Kohama, Y. 1987 Three-dimensional boundary layer transition on a convex- 
concave wall. Paper presented at turbulence meeting , Bangalore , January 1987. 

Peerhossaini, H. fc Wesfreid, J.E. 1988 On the inner structure of streamwise 
Gortler rolls. Inti J . Heat Fluid Flow, 9, 12-18, 

Smith, A.M.O. 1955 On the growth of Taylor- Gortler vortices along highly concave 
walls. Q. Appl. Math., 13, 233-262. 

Wortmann, F.X. 1969 Visualisation of transition. J . Fluid Mech. 38, 473-480. 


26 



FIGURE CAPTIONS 


Figure 1. The different regions beyond the downstream position of neutral stabil- 
ity. 


Figure 2. Eigenfunctions for the weakly nonlinear problem of §4 as defined by 
(4.1)— (4.4) when (A”i)i, = 0.690, tli = 0.372. (a) Re(B\o)', ( b ) iZe(Cio); (c) Re(G\)\ 
( d ) Re(G 2 ); (c) An (/) Re(B 2 o) and (g) Re{Dio). 


Figure 3. Eigenfunctions for the weakly nonlinear problem of §4 as defined by 
(4.1)-(4.4) when (Ki) l = 2.900, (l L = 0.659. (a) IZe(Bio); ( b ) Re{C 10 ); (c) Re(Gi); 
(d) Re(G 2 ); (e) Ai; (/) Re{B 2 o) and ( g ) Re(Dio). 
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